On the effects of solenoidal and compressive turbulence in prestellar cores

We present the results of an ensemble of SPH simulations that follow the evolution of prestellar cores for $0.2\,{\rm Myr}$. All the cores have the same mass, and start with the same radius, density profile, thermal and turbulent energy. Our purpose is to explore the consequences of varying the fraction of turbulent energy, $\delta_\mathrm{sol}$, that is solenoidal, as opposed to compressive; specifically we consider $\delta_\mathrm{sol}=1,\,2/3,\,1/3,\,1/9\;{\rm and}\;0$. For each value of $\delta_\mathrm{sol}$, we follow ten different realisations of the turbulent velocity field, in order also to have a measure of the stochastic variance blurring any systematic trends. With low $\delta_\mathrm{sol}(<\!1/3)$ filament fragmentation dominates and delivers relatively high mass stars. Conversely, with high values of $\delta_\mathrm{sol}(>\!1/3)$ disc fragmentation dominates and delivers relatively low mass stars. There are no discernible systematic trends in the multiplicity statistics obtained with different $\delta_\mathrm{sol}$.Comment: 9 pages. Accepted by MNRA

Filamentary fragmentation in a turbulent medium

We present the results of smoothed particle hydrodynamic simulations investigating the evolution and fragmentation of filaments that are accreting from a turbulent medium. We show that the presence of turbulence, and the resulting inhomogeneities in the accretion flow, play a significant role in the fragmentation process. Filaments which experience a weakly turbulent accretion flow fragment in a two-tier hierarchical fashion, similar to the fragmentation pattern seen in the Orion Integral Shaped Filament. Increasing the energy in the turbulent velocity field results in more sub-structure within the filaments, and one sees a shift from gravity-dominated fragmentation to turbulence-dominated fragmentation. The sub-structure formed in the filaments is elongated and roughly parallel to the longitudinal axis of the filament, similar to the fibres seen in observations of Taurus, and suggests that the fray and fragment scenario is a possible mechanism for the production of fibres. We show that the formation of these fibre-like structures is linked to the vorticity of the velocity field inside the filament and the filament's accretion from an inhomogeneous medium. Moreover, we find that accretion is able to drive and sustain roughly sonic levels of turbulence inside the filaments, but is not able to prevent radial collapse once the filaments become supercritical. However, the supercritical filaments which contain fibre-like structures do not collapse radially, suggesting that fibrous filaments may not necessarily become radially unstable once they reach the critical line-density.Comment: (Accepted for publication in MNRAS

Star Formation triggered by cloud-cloud collisions

We present the results of SPH simulations in which two clouds, each having mass $M_{_{\rm{o}}}\!=\!500\,{\rm M}_{_\odot}$ and radius $R_{_{\rm{o}}}\!=\!2\,{\rm pc}$, collide head-on at relative velocities of $\Delta v_{_{\rm{o}}} =2.4,\;2.8,\;3.2,\;3.6\;{\rm and}\;4.0\,{\rm km}\,{\rm s}^{-1}$. There is a clear trend with increasing $\Delta v_{_{\rm{o}}}$. At low $\Delta v_{_{\rm{o}}}$, star formation starts later, and the shock-compressed layer breaks up into an array of predominantly radial filaments; stars condense out of these filaments and fall, together with residual gas, towards the centre of the layer, to form a single large-$N$ cluster, which then evolves by competitive accretion, producing one or two very massive protostars and a diaspora of ejected (mainly low-mass) protostars; the pattern of filaments is reminiscent of the hub and spokes systems identified recently by observers. At high $\Delta v_{_{\rm{o}}}$, star formation occurs sooner and the shock-compressed layer breaks up into a network of filaments; the pattern of filaments here is more like a spider's web, with several small-$N$ clusters forming independently of one another, in cores at the intersections of filaments, and since each core only spawns a small number of protostars, there are fewer ejections of protostars. As the relative velocity is increased, the {\it mean} protostellar mass increases, but the {\it maximum} protostellar mass and the width of the mass function both decrease. We use a Minimal Spanning Tree to analyse the spatial distributions of protostars formed at different relative velocities.Comment: 10 pages, 11 figure

Star formation triggered by non-head-on cloud-cloud collisions, and clouds with pre-collision sub-structure

In an earlier paper, we used smoothed particle hydrodynamics (SPH) simulations to explore star formation triggered by head-on collisions between uniform-density 500 M clouds, and showed that there is a critical collision velocity, vCRIT. At collision velocities below vCRIT, a hub-and-spoke mode operates and delivers a monolithic cluster with a broad mass function, including massive stars (M 10 M) formed by competitive accretion. At collision velocities above vCRIT, a spider’s-web mode operates and delivers a loose distribution of small sub-clusters with a relatively narrow mass function and no massive stars. Here we show that,if the head-on assumption is relaxed, vCRIT is reduced. However, if the uniform-density assumption is also relaxed, the collision velocity becomes somewhat less critical: a low collision velocity is still needed to produce a global hub-and-spoke system and a monolithic cluster, but, even at high velocities, large cores – capable of supporting competitive accretion and thereby producing massive stars – can be produced. We conclude that cloud–cloud collisions may be a viable mechanism for forming massive stars – and we show that this might even be the major channel for forming massive stars in the Galaxy

First Investigation of the Combined Impact of Ionizing Radiation and Momentum Winds from a Massive Star on a Self-gravitating Core

J. Ngoumou, et al., “First Investigation of the Combined Impact of Ionizing Radiation and Momentum Winds from a Massive Star on a Self-gravitating Core”, The Astrophysical Journal, Vol. 798(1), December 2015. © 2015. The American Astronomical Society.Massive stars shape the surrounding interstellar matter (ISM) by emitting ionizing photons and ejecting material through stellar winds. To study the impact of the momentum from the wind of a massive star on the surrounding neutral or ionized material, we implemented a new HEALPix-based momentum-conserving wind scheme in the smoothed particle hydrodynamics (SPH) code SEREN. A qualitative study of the impact of the feedback from an O7.5-like star on a self-gravitating sphere shows that on its own, the transfer of momentum from a wind onto cold surrounding gas has both a compressing and dispersing effect. It mostly affects gas at low and intermediate densities. When combined with a stellar source's ionizing ultraviolet (UV) radiation, we find the momentum-driven wind to have little direct effect on the gas. We conclude that during a massive star's main sequence, the UV ionizing radiation is the main feedback mechanism shaping and compressing the cold gas. Overall, the wind's effects on the dense gas dynamics and on the triggering of star formation are very modest. The structures formed in the ionization-only simulation and in the combined feedback simulation are remarkably similar. However, in the combined feedback case, different SPH particles end up being compressed. This indicates that the microphysics of gas mixing differ between the two feedback simulations and that the winds can contribute to the localized redistribution and reshuffling of gas.Peer reviewe

Smoothed Particle Hydrodynamics simulations of expanding HII regions. I. Numerical methods and tests

We describe a new algorithm for including the dynamical effects of ionizing radiation in SPH simulations, and we present several examples of how the algorithm can be applied to problems in star formation. We use the HEALPix software to tessellate the sky and to solve the equation of ionization equilibrium along a ray towards each of the resulting tesserae. We exploit the hierarchical nature of HEALPix to make the algorithm adaptive, so that fine angular resolution is invoked only where it is needed, and the computational cost is kept low. We present simulations of (i) the spherically symmetric expansion of an HII region inside a uniform-density, non--self-gravitating cloud; (ii) the spherically symmetric expansion of an HII region inside a uniform-density, self-gravitating cloud; (iii) the expansion of an off-centre HII region inside a uniform-density, non--self-gravitating cloud, resulting in rocket acceleration and dispersal of the cloud; and (iv) radiatively driven compression and ablation of a core overrun by an HII region. The new algorithm provides the means to explore and evaluate the role of ionizing radiation in regulating the efficiency and statistics of star formation.Comment: 12 pages, 16 figures, simulation movies available at http://galaxy.ig.cas.cz/~richard/HIIregion

GANDALF - Graphical Astrophysics code for N-body Dynamics And Lagrangian Fluids

GANDALF is a new hydrodynamics and N-body dynamics code designed for investigating planet formation, star formation and star cluster problems. GANDALF is written in C++, parallelised with both OpenMP and MPI and contains a python library for analysis and visualisation. The code has been written with a fully object-oriented approach to easily allow user-defined implementations of physics modules or other algorithms. The code currently contains implementations of Smoothed Particle Hydrodynamics, Meshless Finite-Volume and collisional N-body schemes, but can easily be adapted to include additional particle schemes. We present in this paper the details of its implementation, results from the test suite, serial and parallel performance results and discuss the planned future development. The code is freely available as an open source project on the code-hosting website github at https://github.com/gandalfcode/gandalf and is available under the GPLv2 license.This research was supported by the DFG cluster of excellence "Origin and Structure of the Universe", DFG Projects 841797-4, 841798-2 (DAH, GPR), the DISCSIM project, grant agreement 341137 funded by the European Research Council under ERC-2013-ADG (GPR, RAB). Some development of the code and simulations have been carried out on the computing facilities of the Computational centre for Particle and Astrophysics (C2PAP) and on the DiRAC Data Analytic system at the University of Cambridge, operated by the University of Cambridge High Performance Computing Service on behalf of the STFC DiRAC HPC Facility (www.dirac.ac.uk); the equipment was funded by BIS National E-infrastructure capital grant (ST/K001590/1), STFC capital grants ST/H008861/1 and ST/H00887X/1, and STFC DiRAC Operations grant ST/K00333X/1

Protoplanetary disc evolution affected by star-disc interactions in young stellar clusters

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society. © 2014 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society.Most stars form in a clustered environment. Therefore, it is important to assess how this environment influences the evolution of protoplanetary discs around young stars. In turn, this affects their ability to produce planets and ultimately life. We present here for the first time 3D smoothed particle hydrodynamics/N-body simulations that include both the hydrodynamical evolution of the discs around their natal stars, as well as the dynamics of the stars themselves. The discs are viscously evolving, accreting mass on to the central star and spreading. We find penetrating encounters to be very destructive for the discs as in previous studies, although the frequency of such encounters is low. We also find, however, that encounter influence the disc radii more strongly than other disc properties such as the disc mass. The disc sizes are set by the competition between viscous spreading and the disruptive effect of encounters. As discs spread, encounters become more and more important. In the regime of rapid spreading, encounters simply truncate the discs, stripping the outer portions. In the opposite regime, we find that the effect of many distant encounters is able to limit the disc size. Finally, we predict from our simulations that disc sizes are limited by encounters at stellar densities exceeding ∼2–3 × 103 pc−2.Peer reviewe